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Joint measurements of non-commuting observables are characterized by unavoidable measurement 
uncertainties that can be described in terms of the error statistics for input states with well-defined 
values for the target observables. However, a complete characterization of measurement errors 
must include the correlations between the errors of the two observables. Here, we show that these 
correlations appear in the experimentally observable measurement statistics obtained by performing 
the joint measurement on maximally entangled pairs. For two-level systems, the results indicate 
that quantum theory requires imaginary correlations between the measurement errors of X and 
Y since these correlations are represented by the operator product XY = iZ in the measurement 
operators. Our analysis thus reveals a directly observable consequence of non-commutativity in the 
statistics of quantum measurements. 
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I. INTRODUCTION 

Recently, advances in experimental possibilities have renewed the interest in the physics of measurement uncer¬ 
tainties This topic of research actually has a long history going all the way back to Heisenberg’s justification 

of the quantum formalism by the definition of uncertainty limits for the simultaneous determination of position and 
momentum. Consequently, most of the recent work has focused on the quantitative evaluation of uncertainties using 
either the original evaluation by variances jTol - !T^| or the more recent concept of entropic uncertainties fl3l - il5l |. Cru¬ 
cially, all of these approaches are based on the statistical limitations imposed by the mathematics of the quantum 
formalism, following Heisenberg’s implicit suggestion that it is impossible to obtain any direct experimental evidence 
of the physics that causes the appearance of uncertainties in the first place. However, the study and characterization 
of entanglement shows that quantum correlations can exceed the local limits of uncertainties (16l-H8i|. It is therefore 
possible to observe otherwise hidden details of the measurement error statistics by using the uncertainty free correla¬ 
tions of entangled states as a reference. In the present paper, we show how the correlations between the errors that 
characterize a joint measurement of two non-commuting observables can be obtained by analyzing the experimental 
results from maximally entangled inputs. Specifically, we show that the reconstruction of error statistics is indepen¬ 
dent of the model used for the measurement process or for the physics of the input state, resulting in an operational 
procedure for the reconstruction of error statistics. For two-level systems, the complete statistics of measurement er¬ 
rors is obtained. Significantly, the error statistics for quantum measurements consistent with the standard formalism 
of measurement theory exceed the bounds for real-valued measurement errors, effectively resulting in a derivation of 
complex-valued error probabilities from the experimentally observed data. 

The starting point of our discussion is a specific form of a joint measurement of the non-commuting observables 
A and B , where the results of the joint measurement correspond to the cP combinations of eigenvalues observed in 
separate projective measurements of the two observables. For this kind of measurement, the error statistics of one of 
the two observables is given by the conditional probabilities that are observed experimentally for eigenstate inputs 
of that observable. Measurement errors can then be described in qualitative terms, based on the notion that the 
outcome is either correct or incorrect. The marginal error probabilities for either A or B are directly observed in 
measurements of the eigenstate inputs, while the correlations between errors need to be reconstructed from inputs 
with known correlations. By representing the measurement in terms of a joint error probability, we are essentially 
comparing the quantum statistics of measurements directly with the analogous classical statistics, similar to the way 
in which the Wigner function describes quantum statistics as an analog of classical phase space distributions. The 
problem we need to solve is that the correlations between non-commuting properties in a general quantum state are 
themselves unknown. For this purpose, we introduce the entangled states as an input that probes the correlations 
between measurement errors of non-commuting observables. Specifically, we make use of the fact that maximal 
entangled states provide precise correlations between the values of all physical properties in the two systems, so that 
an error free measurement results in a well-defined relation between the outcomes in system 1 and the outcomes in 
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system 2 for both A and B. It is then possible to deduce the error statistics using only the assumption that the joint 
probabilities for relations other than the ones observed in precise measurements of either A or B are zero. 

We apply our analysis of measurement statistics to the orthogonal Bloch vector components X and Y of a two 
level system and find that the only missing element in the description of the measurement process is the correlation 
between the errors in X and Y. To obtain this correlation experimentally, we then introduce a maximally entangled 
state of two systems. Since the correlations of X and Y between the two systems are known, the application of joint 
measurements to both systems should ideally produce the same correlations. From the actual results, we can judge 
whether the occurrence of an error in only one of the two measurements changed the correlation or not. Since the 
changes in correlations are obtained for both X and Y in a single measurement, the measurement outcomes reveal 
the correlation between the occurrence of errors in the two local measurements, and the missing element of the error 
statistics can be derived from the experimental results. 

We can use the standard formalism of quantum mechanics to predict the results of joint measurements on maximally 
entangled states. Interestingly, the predicted experimental results violate the limits that apply to any positive valued 
(and hence classical) selection of error probabilities. In fact, the quantum result indicates that the error correlation 
must be given by an imaginary part in the error probabilities. When these complex error statistics are applied to 
reconstruct the quasi-probabilities of an arbitrary input state, the result is the Kirkwood-Dirac distribution that 
is also observed in weak measurements mum. Since this distribution corresponds to the expectation value of an 
ordered product of the projection operators, the imaginary parts are a direct consequence of non-commutativity 
and essentially identify the expectation values of commutation relations with imaginary statistical correlations. The 
experimentally observed correlations between the measurement outcomes for the local measurements thus indicate 
that non-commutativity describes imaginary statistical correlations between real-valued observables. 

To properly appreciate the significance of the results, it is necessary to keep in mind that the physics of measurement 
can only be understood in terms of experimentally observable evidence. The essential insight of the following discussion 
is that entanglement provides us with the necessary tool for the analysis of correlations between non-commuting 
properties. If the claim that non-positive probabilities are a necessary consequence of the experimental evidence 
seems to be somewhat odd, it should be kept in mind that the observation or preparation of a joint reality for the 
non-commuting properties in question is impossible. The following discussion is therefore aimed at connecting the 
experimental evidence with the physical properties in a way that identifies the relevant correlations without requiring 
any unphysical assumptions about the reality of the system. 


II. ERROR STATISTICS OF JOINT MEASUREMENTS 

The established formalism of quantum mechanics represents physical properties by operators in Hilbert space. 
The precise measurement of a physical property A is represented by a projection onto an eigenstate | a), where the 
probability of an outcome a is given by the product trace of the projection operator | a) (a | and the density matrix 
p = | ip) {ip I that described the state | ip) of the system before the measurement. Importantly, a joint assignment 
of outcomes a and b for two physical properties A and B is only possible if the two physical properties have shared 
eigenstates. It is therefore impossible to define a joint measurement of a and b without introducing some form 
of measurement uncertainty, such that the probability P(a, b) of a joint outcome is not given by an intrinsic joint 
probability p{a , b) of the initial state p. Nevertheless, we can design a joint measurement that yields a complete set 
of joint outcomes (a, b) if we accept statistical errors in the joint measurement. 

In standard quantum theory, such a joint measurement is described by a positive operator-valued measure, II a ,b) 
so that the joint probability of the experimental measurement outcome is given by 

P exp .(a, b) = Tr{fl ayb p}. (1) 

Mathematically, this is a bilinear relation between the measurement and the state defined in the Hilbert space of 
the system. However, it is the purpose of the measurement to evaluate the physical properties A and B in terms of 
their precise measurement outcomes a and b that would have been obtained from the original input state p. For this 
purpose, we should express the input state in terms of the set of outcomes (a, b) associated with precise error-free 
measurements. Here, it is interesting to observe that the quantum state p in a d-dimensional Hilbert space is described 
by d 2 independent parameters, corresponding to the number of possible combinations of measurement outcomes a 
and b. Therefore, the only possible form of a bilinear relation of the measurement n aj f, and the state p expressed in 
terms of the outcomes a and b is given by 

Pexp.(a, b) = ^2 -P( a > b \ a 'i b ')p( a 'i b '), 

a' ,b' 


( 2 ) 
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where the representation of the quantum state correspond to a joint probability p(a, b) of the outcomes a and 6, and the 
description of the measurement process corresponds to a conditional probability relating the input combination (a ', b') 
to the measurement result (a, b). Importantly, Eq.© is not based on a hidden variable model and does not require 
the assumption that the combination of a and b in p(a , b ) represent a joint reality of a and b before the measurement 
is performed. The purpose of jointly assigning a and b to the input is merely to obtain a formulation of the joint 
measurement II ai h in terms of the target observables A and B that is equally valid for any input state, including 
eigenstates of either A or B. Effectively, Eq.© should be understood as a general formulation of measurement that 
is analogous to classical statistics, just like a Wigner distribution is analogous to a classical phase space distribution. 
Since our goal is the evaluation of the error statistics P(a 1 b\a r ,b'), we will try to make only minimal assumptions 
about the form of the quasi-probability p(a,b). It should be noted that this is possible because the measurement 
errors can be reconstructed from only a limited selection of input states. In the following, the states we consider for 
this purpose are eigenstates of the observables and maximally entangled states. In both cases, the form of the joint 
probabilities p(a, b ) and p(ai, &i; 02 ^ 2 ) is completely determined by the experimental evidence and a set of reasonable 
constraints that would also be valid in classical probability theories. By using entangled states, it is possible to derive 
the non-classical correlations between errors using only the apparently classical correlations observed in separate 
measurements of either A ot B. 

Much like the definition of quasi-probabilities for quantum states, the formulation of P(a, b\a ', b') is motivated 
by the similarity between the description of classical measurement errors and the errors directly observed when the 
input value of either A or B is known. The procedure developed in the following allows a reconstruction of the 
quasi-probability P(a, b\a', b') from experimental results that correspond to marginal probabilities of the complete 
error statistics. Similar to the reconstruction of a Wigner function from marginal distributions, it is natural that a 
reconstruction of the unobservable conditional probability P(a, b\a !, b') will result in non-positive quasi-probabilities. 
Nevertheless it should be kept in mind that these results are obtained from statistical assumptions about the relations 
between physical properties that appear to be completely consistent with all experimental observations. 

The argument is perhaps easiest to understand if we first look at the marginal distributions of P{a , b\a', b') that are 
obtained by performing the measurement on a known eigenstate of A or B. In these cases, it is possible to identify 
a unique joint probability distribution of a and b, since the eigenstates of A or B assign a probability of zero to any 
outcome other than the one specified by the state. The joint probabilities for the correct outcome are then given by 
the marginal probabilities for the eigenstates of the other property, as given by the squared inner products \(b \ a )| 2 . 
For an eigenstate | a), the joint probabilities p a {o !, b') are all zero for a' ^ a and correspond to the probabilities of b' 
in | a) otherwise, 

Pa(a',b') = 5 a y\(b'\a)\ 2 . (3) 

Likewise, the eigenstates of B are described by joint probabilities of 

p b (a',b') = S bM \(b\a')\ 2 . (4) 

Importantly, we assume that there are no hidden non-classical probabilities for values of a'(b') other than the eigenvalue 
a(b). Specifically, a quasi-probability p a (a',b') for an eigenstate | a) could also be formulated in such a way that the 
marginal for a' a is zero because positive and negative quasi-probabilities cancel. Strictly speaking, the assumption 
that the value of A in an eigenstate | a) does not depend on the value of B is therefore one of the constraints that we 
choose to make a reconstruction of the total measurement statistics possible. However, it should be kept in mind that 
we are trying to find the closest possible analogy between quantum statistics and classical statistics, so the assumption 
of any joint probability other than the ones given in Eqs.© and © would seem to be rather artificial, given that 
classical statistics is sufficient to fully explain the experimental results in that specific case. 

We can now characterize the conditional probabilities P(a, b\a r , b') experimentally by measuring the probabilities of 
the outcomes (a, b) for input states of | a") and | b"). The result is an intuitive description of measurement errors in 
terms of conditional probabilities that relate the correct input values to the actual output values. However, the fact 
that we have to select either eigenstates of A or eigenstates of B for the input means that we do not obtain information 
on the detailed correlations between the measurement errors (a" —>• a) and the measurement errors (6" —>• b). The 
experimental data only provides us with the input marginals of P(a, b\a\ b') given by 

Pexp. («,%") = £P(a,%", &')!<&' I a")I 2 , 

b' 

P exp Xa,b\b") = Y / P^b\a',b' , )\(b" |a')| 2 . (5) 

a' 

Eigenstate inputs therefore permit only an incomplete characterization of the measurement statistics. To characterize 
the complete set of correlated errors described by the measurement operator n a b, it is necessary to use input states 
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with well-defined correlations between the non-commuting properties A and B. This is a non-trivial problem, since 
there is no clear consensus on the description of non-classical correlations by quasi-probabilities, and it is usually 
thought that such descriptions can only be obtained in the form of specific statistical models [H, H3]. Fortunately, 
this problem can be sidestepped by using maximally entangled states, and we will show in section ITVl how correlations 
between the measurement errors can be identified in the experimental statistics obtained from correlated measurement 
of maximally entangled pairs. Before moving on to entanglement, however, it may be useful to take a closer look at 
the specific case of two-level systems, where it is sufficient to characterize the measurement outcome as either correct 
or incorrect. 


III. EVALUATION OF MEASUREMENT ERRORS IN TWO-LEVEL SYSTEMS 

The most simple example of a pair of non-commuting observables is give by the orthogonal spin components X 
and Y of a two-level system, where the eigenvalues are given by ±1. A joint measurement of A' and Y therefore 
has four possible outcomes given by ( x,y ) = (+1, +1), (+1, — 1), (—1, +1), (—1, — 1). As explained above, the joint 
measurement is described by 16 conditional probabilities P{x,y\x' ,y') that relate the input values of {x',y') to the 
output values ( x , y). However, the number of unknowns can be greatly reduced if we assume that the error probabilities 
are symmetric under exchanges of +1 and —1, so that the probability of obtaining a correct outcome does not depend 
on the actual spin value. The measurement statistics is then described by only four error probabilities defined by the 
relation between the input values and the output values. We can define these error probabilities as 

7?(0,0) = P{x,y\x,y) 

77 ( 0 , 1 .) = P(x,-y\x,y) 

77 ( 1 , 0 ) = P(—x, y\x, y) 

= P(-x,-y\x,y), (6) 

where the arguments of the error probabilities 77 (r x ,r y ) define whether an error occurs (ri = 1) or not (?y = 0). 
Experimentally, it is always possible to satisfy this condition by randomly flipping the spin direction in the input, 
so that any differences in the errors for the four different input combinations of x and y average out in the overall 
measurement statistics. The removal of any experimental bias in x or 77 is therefore only a technical problem with a 
sufficiently simple solution. 

Following the procedure outlined in section ITT1 we can now characterize the error probabilities ri(r x ,r y ) from exper¬ 
imental data obtained with eigenstate inputs. For the A-eigenstates, the experimental results are 

P exp . (x, y\x) = ^ ( 77 ( 0 ,0) + 77 ( 0 ,1)) 

P exp .(-x,y\x) = ^ (?7(1,0)+ 77 ( 1 , 1 )) (7) 

and these two results can be summarized by a single resolution parameter, 

14 = ^ \ (Pexp. ( X , y\x) Pexp. ( X, y\x)) 

V 

= 77 ( 0 , 0 )+ 77(0,1)- 77 ( 1 , 0 )- 77 ( 1 , 1 ). (8) 

The resolution parameter V x represents the statistical contrast or visibility between the correct outcome of x and the 
opposite outcome. 14 = 1 denotes a precise measurement and 14=0 describes a completely random measurement 
outcome, uncorrelated with the input. The results obtained from Y-eigenstates can be summarized in the same 
manner, by 


14 — ^ ' (Pexp. (x, y\y) Pexp. (x, y\y)) 

X 

= 77 ( 0 , 0 )- 77 ( 0 , 1 )+ 77 ( 1 , 0 )- 77 ( 1 , 1 ). (9) 

The resolution parameters 14 and V y can be evaluated directly by performing the measurement on eigenstate inputs 
and taking the difference between the probability of the correct measurement result and the probability of an error 
in the measurement of the known input observable. In addition, normalization requires that 


77(0,0) + 77(0, 1) + 77(1, 0) + 77(1,1) = 1- 


( 10 ) 


This leaves us with a single unknown parameter, which can be expressed by 

C = 77 ( 0 ,0) - 77 ( 0 ,1) - 77 ( 1 ,0) + 77 ( 1 , 1 ). 
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( 11 ) 

This expression describes the correlation between errors in X and errors in Y by distinguishing whether the errors 
occur jointly or separately. To evaluate this parameter, it is therefore necessary to know the correlations between X 
and Y in the input state, and this makes it necessary to consider the case of a maximally entangled input. 


IV. CORRELATIONS BETWEEN ENTANGLED PAIRS 


Maximally entangled pairs are states with particularly strong correlation between two quantum systems. Specif¬ 
ically, a precise measurement of one system will project the remote system into an eigenstate corresponding to the 
measurement result obtained locally. If the physical properties of the two systems are properly aligned, the value 
of A in the remote system can be determined by an A-measurement on the local system, and the value of B in the 
remote system is determined from the result of a B-nreasurement on the local system. Maximally entangled pairs are 
therefore characterized by correlations in both A and B , providing possible experimental evidence about the quantum 
correlations between the non-commuting observables A and B. 

To obtain a joint probability for entangled pairs, we make use of the fact that the correlations observed in A and 
in B are completely consistent with classical statistics. As for the case of eigenstates given above, we assume that 
no non-positive probabilities should be introduced if the measurement statistics can also be explained in terms of 
conventional positive valued probabilities. In the case of maximally entangled states, this results in a well-defined 
joint probability for the measurement outcomes of A and B. However, it should be noted that the values obtained 
for A and for B in both systems are not normally the same, so it may be important to consider the specific map 
between the values in system 1 and the values in system 2. In the following, we will represent this map by a tilde 
above the variable in system 1, which indicates that the value represented by the expression is actually that observed 
in system 2 as a consequence of the result in system 1, so that 02 = a\ and b 2 = b±. The joint probabilities of precise 
measurements in the two systems for the maximally entangled state | E) can then be given by 

P{ai,a 2 ) = |(ai,a 2 | E )| 2 = i<5 5l! a 2 

P(b ll b 2 ) = \{biM\E)\ 2 = - d 5 iiM . (12) 

It is also possible to determine the probabilities for a measurement of A in one system and a measurement of B in 
the other, 

P(ai,b 2 ) = |(ai, b 2 \ E )| 2 = ^|(ai | b 2 )| 2 

P(b 1 ,a 2 ) = \(b 1 ,a 2 \E)\ 2 = ±\(a 2 \b 1 )\ 2 . (13) 

Given the marginal probabilities above, it is possible to derive the complete joint probability of the maximally 
entangled state without any further assumptions. Specifically, Eq. m indicates that all joint probabilities are zero 
unless 0,2 = ai and b 2 = 61 , and Ea. fTdl) provides the specific probabilities of all non-zero contributions. Therefore, 
there is only one possible quasi-probability describing the maximally entangled state, 


p E (a>i,bi; a 2l b 2 ) = -5- ai ^ 2 5~ bi b2 \{a 2 \ b 2 )| 2 . 


(14) 


With this quasi-probability, it is possible to describe the correlated statistics of two independent measurements 
performed separately on the two systems. Here, we assume that it is possible to perform exactly the same measurement 
independently on system 1 and on system 2. Experimentally, this may be somewhat difficult, since the measurement 
setups may have different imperfections. Ideally, one would use the same setup twice using a time delay that is 
sufficiently long to reset the measurement system. If this is not possible, it would be important to ensure that both 
measurement setups are characterized by the same measurement statistics, perhaps even by randomly exchanging the 
roles of system 1 and system 2. Importantly, it is the goal of the measurement to characterize the joint measurement 
of A and B on a single system, so it is an experimental requirement of the present method of analysis that the 
measurement we wish to analyze can be applied independently to both systems. The experimental statistics of 
such an independent application of the measurement to system 1 and system 2 are described by the corresponding 
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products of the two error probabilities, resulting in a quadratic function of the conditional probabilities P{a, b\a ', b') 
that describe the error statistics of the joint measurement of A and B. Specifically, 

-Pex P .(«ij bi; a 2 , b 2 \E) = ^2 P(ai,bi\a', b')P(a 2 , b 2 \a', b')-\(a' \ b ')| 2 . (15) 

a', b' 

In general, this is a very different sum from the one that determines the errors for eigenstate inputs. In particular, 
the sum runs over all a 1 and all b' of the input, since the input values of A and B are equally unknown. Nevertheless 
the correlations between system 1 and system 2 provide a clear statistical structure to the output, thus revealing 
something about the correlations between the measurement errors in A and in B in the form of correlations between 
the output values of (ai,&i) and (a 2 , b 2 ). 


V. ERROR CORRELATIONS FOR TWO-LEVEL SYSTEMS 

We can now take a closer look at the case of singlet entanglement between a pair of two-level systems. In that case, 
all of the spin components have opposite values, so that Xi = —x\ and y\ = —y\ and the joint probability of Xi and 
Ui can be written as 


pe {%i i yi > x 2 1 y 2 ) xi,x 2 &—yi,y2 ■ (l®) 

The experimental statistics for the correlated outcomes of joint measurements of system 1 and system 2 is given by 

Pexp.( x i>yii x 2 ,y 2 \E) = i ^2 P(xi,yi\x',y')P{x 2 ,y 2 \ -x',-y'), (17) 

x',v' 

where it is assumed that the error statistics of the two measurements are exactly identical. In an actual experiment, 
this could be accomplished by using the same measurement setup twice, with a sufficient time delay between the two 
measurements. Of course, some precautions should be taken to avoid cross-talk between the measurements, since this 
could induce artificial correlations. In setups where sufficient temporal separations between the measurements is not 
an option, care must be taken that the differences in the measurement procedures are completely random and uncor¬ 
related, so that differences between the measurements merely add to the statistical errors described by P{x 1 , yi\x', y'). 
In this context, it may also be worth noting that experimental imperfections in the entanglement source can be over¬ 
come either by subtracting the statistical background noise associated with non-maximal correlations, or by simply 
attributing the errors in the entanglement source to the measurement process. In any case, it should not be too 
difficult to confirm the main results of the present paper experimentally even when the experimental setup does not 
achieve high levels of quantum coherence. 

The sums in Ea. (117l) can be further simplified by using the error probabilities y(r Xl r v ) introduced in section Hill 
This means that the 16 experimental probabilities can also be summarized in terms of only four output patterns, 
depending on whether the observed correlations between (xi,y±) and {x 2 ,y 2 ) correspond to the original correlations 
or not: 


£(0, 0 ) = Pexp. (x, y; -x, -y) 

#(0,1) = P exp .{x,y;-x,y) 

£(1,0) = P e * p .{x,y;x,-y) 

£(!,!) = -Pexp. {x, y, x, y). ( 18 ) 

Each of the four values E(r x ,r y ) can be determined experimentally by applying the joint measurement of x and 
Y to both systems of the entangled pair. According to Ea. (fI71) . these experimental results are related to the error 
probabilities i)(r x ,r y ) by 

V 2 = ( 77 ( 0 ,0) + 77 ( 0 , 1 ) - 77 ( 1 , 0 ) - 77 ( 1 , l)) 2 = 4 (£7(0, 0 ) + £7(0,1) - £7(1, 0 ) - £7(1, 1 )) (19) 

V 2 = ( 77 ( 0 , 0 ) - 77 ( 0 , 1 ) + 77 ( 1 , 0 ) - 77 ( 1 , l )) 2 = 4 (£7(0, 0 ) - £( 0 , 1 ) + £7(1, 0 ) - £7(1, 1 )) ( 20 ) 

C 2 = ( 77 ( 0 ,0) + 77 ( 0 ,1) - 77 ( 1 ,0) - 77 ( 1 , l )) 2 = 4 (£7(0,0) + £(0,1) - £7(1,0) - £7(1,1)) (21) 

While Ea. (fl9l) and Eq. (EUl) merely reproduce the resolutions obtained from the measurements of eigenstate inputs, 

Eq.m provides direct experimental evidence of the correlation between measurement errors in X and Y. 
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VI. QUANTUM THEORY OF JOINT MEASUREMENTS 

It is important to note that the results for the measurement resolutions V x ,V y , and C obtained from the experimental 
data given by the outcome probabilities E{r x ,r y ) according to Eas. (119112ll) are completely independent of the statistical 
model used to explain the measurement process. In particular, Eg. (1211) represents an operational definition of the 
correlation between measurement errors in X and Y that applies equally well to classical and to quantum models. 
However, quantum theory prevents an independent confirmation of the validity of Eq. m by individual measurements 
of appropriately prepared inputs, since there are no joint eigenstates of X and Y that could be used to define the 
relation between the two non-commuting properties in the input. In quantum mechanics, Ea. (121l) thus represents the 
most fundamental operational definition of correlations between measurement errors in joint measurements of X and 
Y. 

Quantum theory actually makes very precise predictions about the measurement outcomes obtained from joint 
measurements of X and Y . Since the measurement operators are defined by self-adjoint matrices in a two-dimensional 
Hilbert space, they can be expressed as linear combinations of the Pauli matrices X , Y, Z and the identity I. The 
requirement of symmetry between results of +1 and results of —1 further limits the choice to equal coefficients with 
variable signs, so the only possible form of the positive operator-valued measure for a joint measurement of X and Y 
is 

n + i, + i = — (i + V x X + V y Y + V z zJ 
n+r,-! = i (/ + V x X - VyY - V z z) 

n_i )+1 = \{i- v x x + v y Y - v z z ) 

= \(i- ~ V v f + v *z) ( 22 ) 

Here, V x and V y correspond directly to the experimentally observed resolutions for eigenstate inputs as discussed in 
sec. m above. Since the correct outcomes should have probabilities greater than 1/2, 14 and V y should be positive 
numbers between zero and one. On the other hand, V z is not related to the experimental evidence discussed in sec. nm 
so it represents an arbitrary parameter that can be either positive or negative. In addition, the quantum formalism 
requires that the operators are positive, so the sum of the squares of V, has a maximal value of one, corresponding to 
projections onto pure states. 

It is now possible to predict the experimental results that can be obtained from correlated measurements of max¬ 
imally entangled pairs characterized by expectation values of (XiX 2 ) = —1, (Y{Y 2 ) = —1 and {Z\Z 2 ) = —1. The 
probabilities of the outcomes are given by the expectation values of the measurement operators, 

-Pexp.fal, 2/2) = (A xi,yi ®ft* 2 , 2 / 2 ) 

= (l - xix 2 V x - yiy 2 Vy - xix 2 yiy 2 V 2 ) . (23) 

As expected, there are only four different results E(r x ,r y ), where r x = 0 for X\ = —x 2 and r y = 0 for y\ = — y 2 . 
Specifically, the four experimentally observable probabilities are given by 

£(o,0) = (1 + V x + V v ~ V z) 

E( 0,1) = ^ (1 + K 2 - V v + V z) 

E(l,0) = ±(l-V x 2 + V y 2 + V z 2 ) 

E( 1,1) = ^ (1 - V x - V y 2 - V 2 ) . (24) 

This result confirms the identification of 14 and V y with the measurement resolutions of X and Y according to Eq.© 
and Eq.©, respectively. Interestingly, the sensitivity to Z now appears as a correlations between the measurement 
errors in X and in Y. However, this correlation is anomalous, since the sign of C 2 becomes negative: 

C 2 = ( 77 ( 0 , 0) - 77 ( 0 ,1) - 77 ( 1 , 0) + 77(1, l)) 2 
= 4 {E{ 0,0) - E( 1,0) - E{ 0,1) + E( 1,1)) 

= —V 2 . 

z 


( 25 ) 
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The experimentally observed correlations between measurement errors in X and Y therefore correspond to imaginary 
error probabilities defined by 


c = (v(0, 0) - 77 ( 0 , 1) - 77 ( 1 ,0) + 77 ( 1 , 1)) = ±iV z . (26) 

This result clearly describes a qualitative difference between the predictions of quantum theory and the predictions 
of any possible classical model of measurement errors. Specifically, all classical statistical models would require that 
£(0,0) + £7(1,1) is greater than £(0,1) + £(1,0). Indeed, it is obvious that a precise measurement of X and Y 
would have to result in the correct correlations, so that only £(0,0) would obtain a non-zero value. The fact that all 
quantum measurements result in £(0,1) + £(1, 0) > £(0, 0) + £(1,1) is therefore closely linked to the impossibility 
of uncertainty free joint measurements. 


VII. OPERATOR CORRELATIONS AND IDEAL MEASUREMENTS 

In order to understand why quantum theory can produce results that contradict classical expectations in a qualita¬ 
tive way, it is necessary to remember that quantum mechanics does not permit an uncertainty free joint measurement 
of the non-commuting properties X and Y. If experimentally observable relations between non-commuting probabil¬ 
ities violate expectations associated with a hypothetical joint reality, the most likely conclusion is that there is no 
such reality. 

The present analysis allows us to identify the actual measurement errors in a quantum measurement, using a close 
analogy with classical error probabilities. The ideal classical measurement would be characterized by 77 ( 0 , 0) = 1 and 
V x = V y = C = 1. Interestingly, we can construct an operator that represents this ideal classical measurement by 
using V z = ±i in Eg. (1771) . The resulting operators can then be factorized into a product of two projectors, 

riideal = i (j ± 1) (i ± y). (27) 

Quantum theory thus provides a well-defined theoretical form for uncertainty-free measurements, but this form is 
unphysical because it is non-hermitian and would therefore result in complex probabilities for the joint outcomes. 
Specifically, the correlations between X and Y are imaginary and satisfy the operator equation 

(XY)=i(Z). (28) 

Thus, our analysis suggests that (a) the correlation between X and Y is correctly represented by an ordered product 
of the operators, and hence has an imaginary value and (b) error free joint measurements are impossible because 
imaginary error probabilities are needed to convert the imaginary correlations into real-valued probabilities. 

Importantly, our analysis shows that the non-classical error correlations of joint measurements can be observed 
directly in experiments using entangled state inputs, without any prior assumptions about non-connnutative observ¬ 
ables. The correlations between the outcomes of joint measurements performed on entangled pairs thus show that the 
imaginary correlations represented by operator products do have experimentally observable consequences and should 
be taken seriously as valid descriptions of non-classical correlations. As Eg. (1*271) indicates, it is possible to invert the 
relation between experimental probabilities and the quasi-probability p{x, y) to obtain the only possible definition of 
a joint probability for a two level system that is consistent with the experimentally observed statistics. This joint 
probability has the form of a Kirkwood-Dirac distribution HMH, 

P(x,y) = Tr (llideaip) = (x | y)(y \p\x). (29) 

Significantly, this distribution has been derived by analyzing the errors of measurements for specific situations in which 
the form of p(x, y) is not ambiguous and can therefore be formulated without having to chose between different quasi¬ 
probabilities. Our results therefore indicate that the only joint probability that is consistent with the experimentally 
observed correlations between measurement errors is the Kirkwood-Dirac distribution given by Eg. (1291) . 

Consistency requires that the Kirkwood-Dirac distribution satisfies all of the assumptions we made when identifying 
the joint probability distributions of eigenstates and entangled states given in Eqs. (13I4I15I17I) . It is comparatively 
easy to confirm this for the case of the eigenstates. For entangled pairs, the Kirkwood-Dirac distribution is given by 

PE(ai,b 1 ;a 2 ,b 2 ) = (oi,a 2 | bi,b 2 )(bi,b 2 | £)(£ | ai,a 2 ) 

= 1(02 | 62)| 2 ^6. ,6^5! , 02 - 


(30) 
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Here, the delta functions originate from the same inner products of state vectors that also describe the experimentally 
observed statistics given in Eq. m- The Kirkwood-Dirac distribution therefore provides a description of entanglement 
that is particularly close to the classical explanation of the correlations between the properties A and B in the two 
systems. It may be worth noting that this result seems to be consistent with theoretical considerations about the 
conditions that define quasi-probabilities, where it was shown that the Kirkwood-Dirac distribution can be derived 
from a relatively small set of reasonable conditions [HJ . It may well be that the Kirkwood-Dirac distribution provides 
the best expression of the quantum correlations described by the Hilbert space formalism. 

In the specific case of a two level system, the reconstruction of the Kirkwood-Dirac distribution from a joint 
measurement of x and y is particularly simple since it only involves the coefficients Vj- The reconstruction procedure 
essentially involves an inversion of these visibilities, where C = iV z is the imaginary visibility of the correlation. The 
inversion can then be given by 

p^,y) = l(l + ± + ±- + ^jP exp X X ,y) 

If 1 1 1 \ „ 

+ 4 { + Vf~Vy~Cj exp(X ’ ~ y) 

If 1 1 1 \ , 

+ 4 \ 1 ^Vf + Vy~c) Pex P - (-a; ’ y) 

1 / 1 1 1 \ , 

+ 4V 1_ T4 ~Vy + Cj P ^ pi ~ X, ~ y) - (31) 

The coefficients V x , V y and C are directly obtained from the experimental data of entangled states according to 
Eqs. (| 19J) . (|20J) and m, respectively. Since C = iV z is imaginary in two level quantum systems, Ea- lETtl) usually 
results in imaginary parts of the quasi-probability p(x, y ) when applied to the real and positive probabilities of the 
experimental measurement outcomes for x and y. The experimental evidence obtained from entangled state inputs 
therefore indicates that the real and positive probabilities of joint measurement outcomes originate from a combination 
of complex joint probabilities with complex error probabilities, where the imaginary parts of the error probabilities 
are needed to compensate the imaginary correlations that are part of the non-classical statistics of the initial quantum 
states. 

Although it is a widespread opinion that theoretical assumptions are necessary to formulate statistical theories of 
quantum mechanics, the present results suggest that these assumptions can be kept to a very reasonable minimum by 
considering the experimentally observable physics in detail. It is then possible to identify non-classical correlations 
directly in the experimental data, without the need for a specific statistical model. Interestingly, this is also consistent 
with recent progress in the analysis of quantum correlations using weak values, in particular the realization of direct 
measurements of the Kirkwood-Dirac distribution (25l - l27i |. which show that the Kirkwood-Dirac distribution is the 
only quasi-probability that correctly describes the non-classical correlations observed in weak measurements, and the 
observation that optimal cloning maps the unobservable correlations described by the Kirkwood-Dirac distribution 
onto observable correlations between the cloned outputs HE 13 ■ Given all of these results, it may well be possible 
that quantum mechanics can be explained completely by the experimental evidence obtained in properly evaluated 
measurements, and that interpretational difficulties are merely a result of misunderstandings caused by unnecessary 
abstractions in the mathematical formulation. In particular, it has already been shown by one of the present authors 
that the Hilbert space formalism can actually be derived completely from the experimental evidence obtained in weak 
measurements [30(. The present results provide additional confirmation that these correlations are physical, and can 
be obtained from the experimental evidence without any arbitrary theoretical assumptions. 


VIII. CONCLUSIONS 

Even though it is impossible to obtain uncertainty free values of non-commuting observables in a single joint 
measurement, it is possible to describe the statistical errors of joint measurements in terms of conditional probabilities 
relating the error-free results to the actual results. Since these statistical errors should not depend on the input state, 
they can be evaluated by using specific input states whose joint probabilities are known from the experimentally 
accessible statistics. These states are the eigenstates of the observables and maximally entangled states with perfect 
correlations in the two non-commuting properties. For two-level systems, the complete set of error probabilities can 
be obtained by comparing the correlations between the joint outcomes obtained for the two entangled systems with 
the known correlations of the initial entangled state. Since the analysis of the experimental results does not depend on 
any prior assumptions about the measurement statistics, it is possible to compare the predictions of quantum theory 
directly with the corresponding predictions of classical models of measurement errors. 
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Significantly, the experimentally observable correlations predicted by quantum theory are qualitatively different 
from the predictions of classical statistical models, since the correlations between the errors in X and Y correspond 
to imaginary error probabilities. This experimentally observable difference between quantum theory and classical 
statistics is a direct consequence of non-commutativity, since the imaginary probability can be traced to the operator 
product XY = iZ. Uncertainty limited joint measurements of non-commuting observables are therefore sensitive to 
the non-classical correlations described by the non-commutativity of the operators that represent physical properties 
in the quantum formalism. By performing the same type of joint measurement on two maximally entangled systems, it 
is possible to make these non-classical correlations visible in the form of joint probability distributions that cannot be 
explained in terms of positive-valued statistics. The correlations between the outcomes of joint measurements observed 
with entangled input states thus provides direct experimental evidence for the non-classical correlations associated 
with non-commutativity, indicating that non-commutativity corresponds to imaginary correlations between physical 
properties. By taking these experimentally observable features into account, it may be possible to develop a complete 
theory of quantum statistics without any arbitrary assumptions from speculative models of reality. 
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